CLNS 99/1631 
CLEO 99-11 



Resonant structure of r ^ Stttt i^r and r ujtiVj- decays 



CLEO Collaboration 
(February 7, 2008) 

Abstract 

The resonant structure of the four pion final state in the decay r — > ^-mr^v-j- is 
analyzed using 4.27 million r+r^ pairs collected by the CLEO II experiment. 
We search for second class currents in the decay r lotti^t using spin-parity 
analysis and establish an upper limit on the non-vector current contribution. 
The mass and width of the p' resonance are extracted from a fit to the r — > 
utti't spectral function. A partial wave analysis of the resonant structure of 
the r Sn-K^Vr decay is performed; the spectral decomposition of the four 
pion system is dominated by the ujtt and aivr final states. 
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I. INTRODUCTION 



The resonant structure of the decay r — > AttUt- has been the subject of experimental 
studies for the past several decades. The first observations of r — > ujhi't- were reported by the 
ARGUS and the CLEO [| Collaborations in 1987. Since then, the ARGUS Collaboration 
estimated the ptttt branching fractions and the ALEPH Collaboration provided improved 
measurements of branching fractions of the r —>■ AttVt- subchannels. No attempt has been 
made, so far, to determine a complete resonant structure of the four pion final state. 

In the Standard Model, r decays to four pions proceed via vector current producing 
hadrons in the = 1~ spin-parity state. At low energies, the hadronic current is dominated 
by vector p, p' and p" resonance form factors which can be approximated by relativistic 
Breit-Wigner functions. The full formalism is discussed in Refs. The parameters of 

= 1~ resonances are of particular interest for phenomenological tests of the Conserved 
Vector Current (CVC) hypothesis that compare multipion final states in r decays and e~^e~ 
annihilations . So far, there have been only two published attempts [p!0| , pA| to extract the 
mass and width of the p' resonance from analyses of the r -k-k^v^ decay and none from 

In this study, we determine the resonant decomposition of the decay r —>■ utiUt and search 
for second class currents via a spin parity analysis. We fit the r unrur spectral function 
to obtain the relative contributions of the p, p' and p" resonances, as well as the mass and 
width of the p' meson. We assume that in a first approximation the widths of the p, p' 
and p" resonances do not depend on the invariant mass of four pions. We use the extracted 
values of the mass and width of the p' to perform a full unbinned maximum likelihood fit of 
the 27r~7r"'"7r*^ final state. 

To perform this fit, we use several phenomenological models. Besides the un channel, we 
allow for aiTT, crp, /op, non-resonant pmr and non-resonant Stttt*^ channels, as well as for their 
various combinations, and fit for their relative contributions to the Stttt" resonant structure. 
This analysis shows clear dominance of the utt and aiir components. 

With a better knowledge of the four pion subchannel structure, we then perform a fit to 
the r — > ujTTUr spectral function assuming a mass- dependent width for the three p resonances. 

We do not use particle identification to separate pions from kaons. A small component 
containing kaons (e.g., K3n) is included in the r background, but its contribution is negli- 
gible. Generalization to charge conjugated reactions and final states is implied throughout 
this paper. 



II. EXPERIMENT AND DATA SELECTION 

We use data from the reaction e~^e^ t^t^ collected by the CLEO II experiment at 
the Cornell Electron Storage Ring (CESR) at or near the energy of the T{4S). The data 
correspond to a total integrated luminosity of 4.68 fb~^ and contain about 4.27 million r+r~ 
pairs. 



CLEO II is a general-purpose solenoidal magnet detector |T2[. The momenta of the 
charged particles are measured by a 67-layer drift chamber tracking system operating inside 
a 1.5 T superconducting solenoid. Photons and electrons are detected in a 7800-crystal 
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Csl electromagnetic calorimeter. Muons are identified using proportional counters placed at 
various depths in the return iron of the magnet yoke. 

The event selection is designed to use the kinematical properties and the low multiplicity 
of r decays to separate them from events copiously produced in two-photon interactions, the 
process e^e~ — > gg, Bhabha scattering, and muon-pair production. We select events with 
four charged tracks in a l-vs-3 topology by requiring that one of the tracks must be at least 
90° away from all the others. When there is more than one combination satisfying these 
criteria, we minimize the acolincarity angle between the momentum of the single track and 
the total momentum of the 3-prong system to select the most back-to-back combination. 
The single track hemisphere is then used as a tag associated with a one-prong r decay to 
ez7gi/^, ihJi^Ut, ttUt- or pz/^. The other three tracks and at least two well identified photons, 
separated by more than 90° from the tag track, are used to reconstruct the hadronic part of 
the signal decay r Snn^Ur — ^ Sn'yyi'r. 

To ensure well modeled acceptance, momenta of the tracks from signal and tag hemi- 
spheres must exceed 0.025£J6^ and O.ObEi^ respectively, where is the energy of the 
colliding electron or positron beam and typically lies between 5.26 and 5.29 GeV. To sup- 
press decay channels with tt+tt" decays, we require that the impact parameter of 
charged tracks with respect to the beam axis must be less than 5 mm. The momentum of a 
tt" candidate is reconstructed from the energies and directions of two photons, constrained 
to the mass of the 7r°. The two photon candidates must produce showers in the barrel part 
of the crystal calorimeter (i.e., |cos^-y| < 0.71, where 9^ is an angle between the photon 
and e+e" beam axis) and be separated by more than 20° from the closest charged track. 
Energy deposition in the calorimeter is required to have a photon-like lateral profile and be 
greater than 75 MeV for the photon with the smaller energy and greater than 120 MeV for 
the photon with the higher energy. We consider only those two-photon combinations that 
are within 2.5 standard deviations from the nominal 7r° mass. To increase statistics for the 
rho tag, we also use showers in the tag hemisphere, detected both in the barrel and endcaps, 
with energy deposition greater than 30 MeV. 

For electron identification, we require that the momentum of the track be greater than 
0.5 GeV/c and that the energy deposited in the calorimeter be consistent with the momentum 
of the electron candidate measured by the drift chamber: \{E/p—1.0)\ < 0.1. We also require 
that the rate of energy loss due to ionization in the drift chamber must be greater than the 
expected rate minus two standard deviations. A track is identified as a muon if it has 
correlated drift and muon chamber hits and traverses at least three absorption lengths of the 
material, when the momentum of the track is below 2 GeV/c, and at least five absorption 
lengths, when the momentum of the track is above 2 GeV/c. To minimize the probability 
of pion-lepton misidentification for both electron and muon tags, we do not allow additional 
showers in the calorimeter with energy deposition greater than 120 MeV that are unmatched 
to any of the charged tracks. A tag track is identified as a pion if it was not identified as 
a lepton and if the invariant mass of the tag hemisphere is less than 0.5 GeV/c^, where 
the invariant mass of the tag hemisphere is defined as the invariant mass of the tag track, 
assuming the pion mass, and all neutral showers in the tag hemisphere. A tag track is 
identified as a rho if it was not identified as a lepton, if the invariant mass of the two most 
energetic showers in the tag hemisphere is consistent with the 7r° mass within 2.5 standard 
deviations and, finally, if the effective mass of the tag track, assuming the pion mass, and 
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of the 71^ candidate lies between 0.5 and 1.2 GeV/c^. To reduce combinatorial and hadronic 
background for the rho tag, we do not allow any showers with energy deposition greater than 
120 (350) MeV in the tag hemisphere if the center of the shower is more (less) than 30 cm 
from the tag track projection onto the calorimeter. 

A typical event, in which the decay r Snn^u^ is tagged by a one-prong r decay, does 
not deposit any significant extra energy in the calorimeter, other than that associated with 
either the 7r° — > 77 decay or with interactions between charged hadrons and Csl crystals. In 
most cases, additional unmatched energy deposited in the calorimeter is a signature of various 
background events having one or more additional vr^'s. To suppress these backgrounds, we 
use the value of the unmatched energy to veto such events. We reject events if there is at least 
one unmatched energy cluster in the calorimeter, which is not used for tt" reconstruction, 
which is more (less) than 30 cm away from the closest projection of the charged track onto 
the calorimeter and has energy deposition more than 75 (350) MeV. 

Further background reduction is based on the fact that, in a r decay, there is always at 
least one neutrino carrying away undetected energy and momentum. Thus, a typical r event 
has a large missing transverse momentum, , whereas a typical two photon event has a small 
missing transverse momentum and small visible energy, Etot ■ The transverse momentum and 
visible energy are measured using the four reconstructed charged tracks, assuming the pion 
masses for the three charged tracks in the signal hemisphere, and one vr" (two vr^'s for the 
rho tag). To remove two-photon background, we use the following selection criteria ||13|| : 
Pt > 0.l{2Ei,m — Etot), Pt > 0.075Ehm and Etot > 0-6Ebm- To reduce background produced 
by Dalitz decays vr" — > e"^e~7 or by 7 conversions in the detector, we reject any event that 
has a well identified electron among the three charged tracks in the signal hemisphere. To 
reduce combinatorial qq background, we impose a requirement on the undetected neutrino 
mass. The neutrino mass can be calculated as: 

= (P, - P,)2 = + g2 + 2|P,||P,| cos^,, - 2Pfe^P, . (1) 

Because we cannot measure cos 9rq, we assume it to be 1 and require that the resulting 
maximum value of the square of the undetected neutrino mass be positive: M^^^^ > 0. In 
Eqn. (0) the value of M^- is taken as 1.777 GeV/c^, the momentum of the r is expressed as 

\Pr\ = \J El^ — M^, the symbols q and Eg represent the invariant mass and energy of the 
four pion final state, and the energy of the r is set equal to the beam energy Ebm- 

The distribution of the Snir^ effective mass, before the cut on the r neutrino mass is 
applied, is plotted in Fig. |^. In the vicinity of the r mass the production of Snir'^ events 
is strongly suppressed by phase space and the signal-to-background ratio becomes small. 
Hence, we accept only those events that have the effective mass of the 37r7r° final state lying 
in the range from 0.9 to 1.7 GeV/c^. 

After all signal selection criteria are applied, we find 25, 374 data events with an estimated 
background of 577 ± 17 events coming from e~^e~ —>■ qq production and 669 ± 13 events from 
other r decays. The background estimates are based on Monte Carlo. The latter background 
source includes contributions from r~ —>■ 27r~7r"'"27r°z/T- and n^n^Ksi'T ~^ 2T^~Ti^Ti^Vr 

decays whose yields are estimated as 1.6% and 0.7%, respectively. The r background also 
includes small Ktitttt^ and KKnir^ components. Their contributions, calculated using the 
values [g of B{t~ ^ R-n+n-n^Ur) = (0.08 ± 0.04)% and B{t- R- K+n-TT^Ur) = 
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(0.069 ±0.030)% for the branching fractions, are found to be neghgible, as well as background 
contamination from other r decays. The combined efficiency for all tags is estimated to be 
6.75%. The measured branching fraction is B{t Snir^Ur) = (4.19±0.10)%, which does not 
include a conservatively estimated 5% systematic error associated with the data analysis. 
This is is in good agreement with the Particle Data Group value of (4.35 ± 0.10)%. 
Contamination due to 27 physics is estimated to be completely negligible. 



III. THE SPECTRAL FUNCTIONS IN r ^ S-Kn^Ur AND r ^ uiriyr DECAYS 

The shape of the 37r7r° effective mass in r decays is distorted by phase space and helicity 
factors. It is useful to examine the spectral function independent of such effects specific to 
r decays. 



The spectral function for r Snn Ur is defined |16 



V^'''' (a) = ^ ' ^ = f2l 

dq GIV^, q{M^ - q^Ym + 2?^) ' ^ ^ 

where q = M(37r7r°), Kd is an element of the CKM matrix, M^- = 1.777 GeV/c^ and Gi;' is 
the Fermi constant. The differential partial width can be represented as: 



(3) 



dq N dq 1927r3 I3{t ez7,zy^) 

Thus, the spectral function is given by: 

3^,0 ^ ^rfiV(g) 1 g(r^37r7rV) 

N dq q{M^-q^Y{M^ + 2q^) B{t -^eUeVr) I'IttV^/ ^' 

To reduce uncertainty caused by the finite bin width, the factor q{M'^ — q^Y{M'^ + 2g^) in 
the denominator is averaged over the bin width. 

In full analogy with V^'^'^° {q), we can define a spectral function for the r ujitUt- decay. 
To extract the component from the Sirn^ final state, we use a technique similar to the usual 
sideband subtraction. For each bin of g = M(37r7r'^), we select events in the u signal region 
(0.76 GeV/c^ < M(7r-7r+7r°) < 0.81 GeV/c^) and fit the combined sidebands (0.83 GeV/c^ < 
M(7r-7r+7rO) < 0.90 GeV/c^) and (0.60 GeV/c^ < M(7r-7r+7rO) < 0.74 GeV/c^) to a second 
order polynomial. The fit parameters are then used to estimate the background under the u 
peak. To estimate the statistical uncertainty of the non-cj component in the signal region, 
we propagate errors using a full covariance matrix. In addition to the statistical error for 
each bin, there is an uncertainty due the sideband definition. To estimate the size of this 
error for each bin, we vary the sideband definition and change the order of the polynomial 
from 2 to 1. 

In Fig. 1^ the shape of the M{ujtt) distribution, which includes both u signal and combi- 
natorial background, is shown with solid squares. The non-cjvr contribution is shown with 
dashed and dotted lines for different sideband definitions, and its average value is shown 
with open squares. For every bin, the error is calculated as a sum in quadrature of both 
statistical and systematic uncertainties. 



7 




0.7 1,2 1,7 2,2 2,7 

M(7t'7T*7tV) [GeV/c'] 

FIG. 1. The distribution of the invariant mass of Stttt^ is plotted for the data (crosses), total 
background (dashed line) and qq background (hatched area), before the cut on the r neutrino mass. 
Further analysis is limited to events between 0.9 and 1.7 GeV/c^. 
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Beginning with the distributions shown in Figs. |^ and 0, we subtract backgrounds, and 
apply efficiency corrections that were determined using a detailed Monte Carlo simulation 
of the r decay process and the detector response. We then employ Eqn. to extract the 
spectral functions V^'^'" (q) and V'^^lq), as well as the non-tuvr contribution y"-"^-'^'^ (^q^j = 
ySmr _ y^-^ (^q^ All three spectral functions are shown in Fig. ^. 



IV. MODEL FITS TO THE r ^ ujirur SPECTRAL FUNCTION 

The spectral function of the component can be expressed in terms of the weak 
form factor F{q'^) as: 



V'^'^iq) 



12tx 



^'^""Fiq) 



(5) 



Here, gp^^.,, is the p utt strong decay coupling, 7^ is a measure of the weak couphng of the 
p resonance to the weak charged current, and 

PM = ^to--(M^-M,)W-(M„ + M,P) 

is the momentum of the lo in the cutt rest frame. 

The choice of the form factor F{q) is largely uncertain. In the region of ~ the 
form factor F{q) is expected to be dominated by low energy resonances. In this analysis 
we use a form for F{q) that allows contributions from the p resonance and its two radial 
excitations p' and p": 

F{q') = BWpiq^) + A, ■ BW^iq^) + ■ BWp.{q') . (7) 

Here, 



BWpiq' 



Ml 



M2 - g2 _ ,qY^ 

is a relativistic Breit-Wigner amplitude normalized to unity at = 0, and 

Ai = ^B^JJL (9) 

are the ratios of the coupling constants for the different p resonances. 

To model the Breit-Wigner shapes (||), we ffist approximate Fp, Fp/ and Fp// by constants 
and do not include any momentum dependence. The parameters of the p(770) resonance 
have been throroughly studied and are known with great precision |T3[. In contrast, the 
parameters of the higher radial excitations p' and p" are poorly known and can be used as 
fit variables. Unfortunately, events above 1.6 GeV/c^ are suppressed by phase space and 
provide little sensitivity for the heavy p", the mass of which is expected to be around 1.7 
GeV/c^. In this study we only fit for the mass and width of p' and keep widths and masses 
of p and p" fixed at their central PDG values |T^: 
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1 500 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




0.90 1,10 1,30 1,50 1,70 

M(a)7T) [GeV/c'] 

FIG. 2. Distributions of the uj-k effective mass for events from the u band (sohd squares) and 
the sidebands (open squares). The u band distribution includes combinatorial background as well. 
Dotted and dashed lines indicate variations in the shape of the combinatorial background due to 
variations in sideband definition. The estimated average background shape, with errors, is shown 
as open squares. 
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0.16 



0.12 - 



0.08 - 
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0.04 - 
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□ con 

O non-con 




0,00 
0,90 



J I I I I I I I I I L 



1.10 1.30 1.50 

q(37T7T^) [GeV/c'] 



1.70 



FIG. 3. Spectral functions calculated for: Stttt" (open squares), lott (solid squares) and non-wTr 
(open circles) components 
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Mp = 770 MeV/c2 , = 151 MeV/c^ ; 
Mpu = 1700 MeV/c^ , r^,, = 235 MeV/c^ . 

The Conserved Vector Current (CVC) hypothesis predicts that the weak form factor and 
the couphng constants are equal to the corresponding electromagnetic quantities. The weak 
coupling ^ — i> uir^ in r decays is expected to be the same as the electromagnetic 
coupling 7* LUTT^ in the reactions e"'"e~ — > 7*(g-y = y/s) — >• cutt" and u — >• 7(g^ = 0)7r°. All 
these couphngs can be presented as gpunr/lp, where 



^rramp ^^^^ 



\ 3r(p ^ e+e-) 



7p 

Thus, one can use the Vector Dominance Model |T8| to cross-check the values of Ai and A2. 
In the VDM, the coupling constants are energy independent, and therefore the set of fitted 
coupling constants should reproduce the experimentally measured value of: 



r(a; ^ 7r°7) = ^aP^ 



^[1 + ^1 + ^2] 

7p 



(0.72 ± 0.04) MeV/c^ (11) 



where P^ is the photon momentum in the u center-of-mass system. 
To fit the data, we use the following spectral models: 

• Model 1: F(g2) oc \BWp^jjo) + A^BWp, + Aa^W^dyoo) P ■ 

• Model 2: F(g2) oc |fiWp(77o) + AiBWp>\'^ . 

• Model 3: F{q^) oc |5iy,(77o) + ^25W,(i7oo) P • 

• Model 4: F{q^) oc |fiWp(770)P • 

The amplitudes Ai and A2 are assumed to be real. We use minimization to fit various 
resonant models to the ujtx spectral function shown in Fig. ^. To obtain the ujti spectral 
function shown in Fig. ^, we divide the uott spectral function shown in Fig. |^ by a value 
of the branching fraction B{uj — » vr+vr"??'') = 0.888 since the four pion final state does 
not include contributions from any other decays of the uj. We verified that the fit results 
are independent of the chosen initial values by repeating the minimization procedure with 
different sets of initial values. 

The following parameters of the p' are obtained from fits to the data with Models 1 and 

2: 

Model 1: Mp, = (1.52 ± 0.01) GeV/c^, Tp, = (0.38 ± 0.04) GeV/c^, 
Model 2: Mp, = (1.53 ± 0.01) GeV/c^, Tp, = (0.43 ± 0.03) GeV/c^. 

Using these results, we calculate weighted averages of the p' parameters: Mp, = (1.523 ib 
0.010) GeV/c^ Tp, = (0.400 ± 0.035) GeV/c^. 

Models 3 and 4 do not provide an acceptable description of the data. Fits to our data 
corresponding to various combinations of the p resonances are displayed in Fig. ^, and the 
results are summarized in Table |. We conclude that the presence of a p' contribution is 
necessary to achieve an acceptable description of the data. Model 2 is favored, but Model 1 
cannot be excluded. 
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FIG. 4. Fits of the ujtt spectral functions with various models. 
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TABLE I. 

Fit results for various combinations of the p,p',p" resonances to the r uji^Vr spectral 
function. Errors are statistical only and are strongly correlated. Errors on the total width 
r(u; 7r°7) are obtained by propagating the entire covariance matrix. 



Model # 


(MeVVc)"i/2 


Ai 


As 


x'/d.o.f. 


MeV 


1 


16.4 ±0.6 


-0.20 ±0.04 


-0.017 ±0.016 


16.1/11 


0.87 ±0.15 


2 


16.1 ±0.6 


-0.24 ±0.02 




16.6/12 


0.78 ±0.10 


3 


22.7 ±0.2 




-0.067 ±0.017 


342/14 


2.3 ± 0.05 


4 


27.2 ±0.2 






940/15 


3.9 ± 0.03 




0.72 ± 0.04 



V. INCLUSION OF MASS-DEPENDENT WIDTHS INTO THE FIT OF THE 

r ^ ujt:v^ SPECTRAL FUNCTION 

We now perform a fit to the r — *■ ujt^Vt- spectrum assuming a mass-dependent width for 
the three p resonances. For this purpose, we use the results of Section [VII Ej We choose the 
model that includes contributions from ujti, pTcn, and non-resonant Stttt^ channels and use the 
extracted values of i?pO^-^o, -Rp-^-^+, Rp+^-^- and -R^^r to perform a numerical integration 
in momentum space in order to obtain a spectral decomposition of the r — > 3'7r'7r''i/^ channel. 
The results of the integration are shown in Table ||. Here and below, we use the value of 
the mass of the p' extracted from the fits to the cutt spectral function as described in the 
previous Section. 



TABLE II. 

Contributions to the spectral functions of the ujtt and ptttt subchannels of r Stttt'^z/^. 







pvrvr 


non-resonant 


7% 


76% 


p(770) 


36% 


0.2% 


/)'(1520) 


12% 


5% 


constructive interference 


47% 


20% 



With a better understanding of the c^tt and pim spectral functions, we perform a fit 
to the r uj-KV-r spectrum assuming mass- dependent widths for the three p resonances. 
Instead of constant-width Breit-Wigner shapes (^, we use mass-dependent total widths 
rj(g^); i = p, p', p". To simulate the mass-dependent widths, we follow the approach adopted 
in Ref. §. 

Since p ^ tttt is the dominant p decay mode, we take into account only this channel for 
the calculation of Tp^q"^): 
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r 



po- 



p^(g^) Y 



(12) 



where F 



pO 



151 MeV and M„ 



770 MeV are the central width and mass of the p resonance 
fixed at their PDG values [15], q is the hadronic four-momentum, PTt{q^) - 



Iq^/A - Ml is 

the pion momentum in the p rest frame, and p^^lM^) is the same momentum evaluated at 
the pole mass of the p resonance. 

For the total width of the p' resonance, we take into account only contributions from 
p' —>■ 717T , p' ujTi and p' —>■ pTCTT modes: 



rp'(?' 



F„. 



p'O 



B{p' 



TTTT 



B{p'^ 
+ B{p' 




pmr) 



(13) 



where Fp/Q and Mp/ are the central width and mass of the p' resonance that are allowed to 
float in the fit, B{p' nn), B{p' un) and B{p' —>■ pTrvr) are relative branching fractions 
normalized with their sum is equal to one, and Puiq^) is the momentum of the uj in the ujtt 
rest frame (|^). To extract the shape of -Fp7r7r(Q'^), we generate a r ptttcUt sample with 
overall phase-space and ujtt contributions eliminated. The form factor is defined as: 



IdN 



(M2 + 2g2)(M2 -g2)2^ 



(14) 



where dN/dq is the number of Monte-Carlo events in a given bin of the invariant mass q. To 
obtain the shape parameters, we fit the form factor to a second order polynomial. 

The choice of the relative branching fractions of p' is largely uncertain The ratio 
B{p' —>■ pTiTi)/B{p' uJTc) can be extracted from our data. The relative contributions of the 
p' resonance into the ujtt and pnir spectral functions, shown in Table ||, are 12% and 5% 
respectively. Hence, the ratio of the branching fractions is estimated as: 

B{p' — » pvrvr) 0.05 B{t /otttz/t-) 



B{p' 



LUTT} 



0.12 B{t 



where the ratio B{t — > p7r7rz/^)/i3(r utiUt-) is given by the results of our multidimensional 
fit described in Section |VII E| . We obtain B{p' pTni)/B{p' ujtt) = 0.61 and vary the 
value of B{p' — > tttt) in the fit to make sure that the outcome depends weakly on the input 
value. For the final result shown in Fig. H, we use the constraint [l^ B{p' 7!'Tr)/B{p' 
con) ~ 0.32. 

To simulate the mass-dependent width of the p", we account only for the p" prni 
decay, since PDG [|T^ suggests that this is the dominant mode. The mass-dependent width 
of the p" is calculated similarly to the p' pinT contribution: 



r,"{q' 



p"0 



Fp^M},) 



(15) 
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where Fp^o = 1700 MeV and Mpn = 235 MeV are the central width and mass of the p" 
resonance fixed at their PDG values |]15|, and Fp^^T.^q'^) has the same shape as for Tp'^q"^). 



Assuming this model, we fit the r — >■ utii't- spectral function to the data. We make two 
fits: 1) the contribution of the p" resonance is set to zero, and 2) the contribution of the p" 
is allowed to fioat in the fit. In both situations the extracted mass of the p' is close to that 
of the p". Hence, we cannot distinguish between the two resonances, and we assume that 
the shape of the spectral function is affected by the p' resonance only. The extracted mass 
and width of the p' are 1.67 ± 0.03 GeV and 0.8 ± 0.1 GeV, respectively. The fit is shown in 
Fig. |. 

We also fit for the p' and p" contributions fixing the masses and widths of p' and p" to 
their PDG values. The obtained fit does not reproduce the shape of the measured spectral 
function, which implies that the given combination of p' and p" parameters does not provide 
an adequate description of the data. 

Finally, we attempt to fit for the relative phases of Ai and A2, i.e., of p' and p" con- 
tributions, respectively. Our study shows that, varying the phases within ±7r/4, we do not 
introduce any significant changes in the shape of the spectral function. Hence, we can argue 
neither in favor of non-zero phases, nor in favor of phaseless models, and we set the phases 
to zero for simplicity. 

We also fit the r —>■ uJirUr spectral function, assuming another parameterization of 
the total width of the p' resonance. Now the p' —>■ pmr contribution is replaced by 
that of p' ai-n because Section |V11E| suggests that the four pion spectrum is domi- 
nated by the ujtx and aivr channels. The last term in Eqn. (pTSf ) is therefore replaced with 
B{p' — > aiTT){Mp' /q){paj^{q'^) /paj^{Mp,)) , where B{p' aivr) is estimated using the results of 
Section [VII E| . The extracted values of the mass and width of the p' are similar to those 
quoted in this Section, i.e., significantly larger than the values obtained with the mass- 
independent fit. 

Because fits assuming mass-dependent widths rely heavily on the choice of the model for 
the width parameterization, we use the mass and width of p' obtained in the fits with the 
mass-independent width (Section |rv|) . The analysis described in this Section is presented 
mostly as a cross-check. 



VI. SEARCH FOR SECOND CLASS CURRENTS IN r ^ ojiri^r DECAY 

The decay r —>■ ujuvr is expected to proceed through the hadronic vector current mediated 
by p, p', p" and higher excitations. If, however, G parity conservation is broken due to second 
class currents |T^, the decay can proceed through a hadronic axial- vector current mediated, 
e.g., by the 6i(1235) resonance. The difference in spin-parity assignments for each of these 
states is reflected in different polarizations of the uj spin and hence in different expected 
angular distributions of cosx- The angle x is deflned as the angle between the normal to 
the uj decay plane and the direction of the fourth pion measured in the uj rest frame, and 
L is the orbital angular momentum of the uon system. The expected forms |T^] of the cos x 



distribution are listed in Table HI 



We form the distribution of cosx for the r ujnUr events in our data. To eliminate 
combinatorial and non-cj background, we perform a sideband subtraction for every bin, in 
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q(con), Gey/c^ 

FIG. 5. Fit of the lott spectral function with mass-dependent width. 
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TABLE III. 

Expected shapes for different spin-parity assignments in the decay r ujitUt-. 





L 


F(cosx) 


l- 


1 


1 — cos^ X 


1+ 





1 


1+ 


2 


1 + 3 cos^ X 


0^ 


1 


cos^ X 



the same way as described in Section |T|. Using Monte Carlo estimates, we subtract small 
uj contributions from r and qq background events, correct for efficiency, and fit the resulting 
distribution, shown in Fig. with the following function: 

F= (l-e)-F^ + e-F2 . (16) 

Here, F^(cosx) = 1 — cos^x is the shape of the dominant vector contribution and 
F^(cosx) = 1 represents the shape which gives the most conservative estimate of the non- 
standard contribution to the r uj-KVr decay. The value of e obtained from the fit is 
consistent with zero within errors: e = (0.08±2.00) x 10~^. After integration, it can be trans- 
lated into an upper limit on the ratio A^'^'^(non- vector current)/A^'^'^(vector current) < 5.4% 
at 90% CL and < 6.4% at 95% CL, where A^'^'^ represents the number of events gen- 
erated through the corresponding current. An analogous study Q by ALEPH yielded 
A^'^'^(non- vector current )/A^'^'^ (vector current) < 8.6% at 95% CL. The solid line in Fig. ^ 
illustrates the agreement of the (1 — cos^x) description and the data. 



VII. FULL UNBINNED LIKELIHOOD FIT TO THE r ^ Svrvr^z^^ RESONANT 

STRUCTURE 

A general method of obtaining a resonant decomposition of a hadronic final state is a 
partial wave analysis of the data via an unbinned maximum likelihood fit. For any data set 
with limited statistics, the analysis is model-dependent since there is a necessary choice of 
amplitudes to be considered. In this study, we consider four models with contributions from 
cuTT, aiTT, ap, /op, pTTTT, and non-resonant Stttt" channels. 

Results depend on the model chosen for the hadronic current parameterization. A general 
description and some details of the formalism are given below. The amplitude description is 
given in the Appendix. 



A. Definition of the fit function 



The differential width of the r decay into the Stttt*^ final state is described by a factorizable 
sequence of decays: t ^ h + Vr followed by /i — > 47r, and can be expressed as: 
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T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



r(ii=i) 




FIG. 6. The distribution of cosx for r toiriyr events after background subtraction and 
efficiency correction is shown as filled squares with error bars. The solid line represents a fit to the 
vector current model = 1^, the dashed line shows the expected shape for = 12", and the 
dotted line shows the = Iq" model. 
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dT{T 3nn^Ur) = ——\M\^dUFS{T Urh)dUFS{h Stttt") , (17) 



where 

2 

is the square of the corresponding matrix element, dLIPS is an element of the phase space 
volume, and L^'^ and H^^ are lepton and hadronic tensors, respectively. 

The lepton tensor L^^ cannot be calculated without information on the neutrino direction. 
To resolve this problem, we use a general form of the matrix element proposed by Mirkes 
and Kiihn which is explicitly averaged over the neutrino direction. In this approach the 
averaged matrix element is expressed as a sum of 16 terms: 

16 

Li^, = 2{Ml-q^)Y,Tm , (19) 

i=l 

where the weights Wi can be calculated from components of the complex hadronic current 
J^, and Li are kinematic parameters which depend only on the four-momenta of the four 
pions in the Stttt^ final state. 

To calculate J^, we use the KORALB Monte Carlo program with the TAUOLA decay pack- 
age The reconstructed four-momenta of the 2ti~ti^ti^ system are boosted into the rest 
frame of the hadronic system which is denoted by the superscript "cm" . They are then pro- 
jected onto the following coordinate system: the x-axis is along the direction of {P^'T + P^JP) , 
the z-axis is perpendicular to the plane defined by the vectors P'^^ and P^!", and the ?/-axis 
is perpendicular to the x — z plane. We found that the calculated weights reveal large un- 
physical fiuctuations unless we require cosip > —0.97, where ip is the angle between the 
direction of the laboratory and the r, as seen from the hadronic rest frame [^]. 

To check whether the calculation of the L^^^H^l, expression is correct, we choose a param- 
eterization of the hadronic current with form factors and coupling constants that provide 
a good description of the data. We generate a large sample of r ^ Stttt^z/^ events accord- 
ing to the phase space and for each event we calculate the corresponding averaged weight: 
L^^^H^y. As a cross check, we use the same parameterization of the hadronic current to 
generate another sample of unit weight events using the standard full TAUOLA simulation. 
The full T Smr^Ur simulation invokes explicit generation of neutrino momenta and subse- 
quent calculation of the precise matrix element L^^H^^. In the limit of large statistics, both 
methods should give the same results. The shapes of submass and total hadronic energy 
distributions generated with both methods agree, though the distributions obtained with 
the averaged matrix element reveal a higher level of fiuctuations. This is due to the fact that 
the averaging over the neutrino direction introduces an additional smearing which increases 
errors. 



B. Likelihood function 

The matrix element averaged over the neutrino direction can be used as a probability 
density function (p.d.f.) to fit the data of the r — Stttt^z/^ decays. The fitting technique is 
a straightforward extension of that used to fit three-body Dalitz plots. 
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We take advantage of the following convenient notation |T3[. The matrix element 
A4{d,f3\9) is described by a model parameter vector 9 which consists of resonant ampli- 
tudes and phases. The matrix element is a function of observable variables a (e.g., hadronic 
four- momenta) and non-observable variables /3 (e.g., neutrino four- momentum) . The matrix 
element has to be integrated over the non-observable variables: 

fs{a\9) = I IMiaJm'dp , (20) 

where fs{d\9) is the probability to observe a final state with a parameter vector a, given the 
model parameters 9. To compare directly with the data distribution, fg must be multiplied 
by the phase space function 0(a) and by the detector acceptance ?](«): 



frnm = r/(a)0(a)/(a|^"') . (21) 

Since we are trying to determine parameters of the matrix element and have all observables 
fixed by data, we follow the usual Bayesian practice and exchange the argument and the 
parameter in the expression for fg^^, thus interpreting it as a probability distribution for 9, 
given the event observables a. 

Since the parameters of the p.d.f are varied in the fit, the normalization of the matrix 
element is no longer constant and must be recalculated on each step of the fit: 

^ ^ ma) ^ md) 

N{9) Jri{a)(p{d)fs{9\a)dd ' 

Due to the presence of background, we must also include a corresponding background p.d.f. 
Fb, and define the total p.d.f. as: 



Rs/bFs + Fb 
Rs/B + 1 



(23) 



where Rs/b is a signal-to-background ratio. The likelihood function for selected events is 
defined as: 

£0) = n F'^^^ = n ^f^^^vici)m . (24) 

events events ^S/ B ~T~ ^ 

To find the maximum of the likelihood, we minimize —2Ln£, where 

^s/bFs 



^S/B-Ts + i'l 
j=l Rs/B + - j=i 



Ln£{9) = Y: Ln ^ ' ^ + ^ Ln(r^(a)0(a)) . (25) 



The term J2 Ln{ri{d)(j){d)) is constant for the purpose of minimizing F and is therefore 
neglected. Thus, we never need to explicitly evaluate r]{d) and 0(a): the amplitude normal- 
ization A^(^) obtained by Monte Carlo techniques takes them into account automatically. 
When parameter terms can be factored out, the computation of the normalization N{9) can 
be speeded up: 
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N{9) = / r]{a)4>{a)f{9\a)da= / r]{d)4>{a 



a 



da 



Y.C,{e) L{a)<f>{a)f,{a)d 



a 



(26) 



where the index j runs over all possible combinations of the fit parameters 6. During the 
fit, the integrals in the above expression are calculated only once and then multiplied by the 
coefficients Cj which depend on the parameters 6 only. 

C. Signal p.d.f. 



As shown in Section [VII B| , all constant factors and functions of observable variables a 
that multiply the signal p.d.f. (^0]) do not change the position of a maximum of the likelihood 
function (|^) in the parameter space 9. Neglecting such contributions, the signal p.d.f ( pOD 
can be expressed as: 



fs = Lf^'H, 
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The hadronic current for the decay r — >■ 27r tt+tt'^z/t- is parameterized as follows: 



(27) 



(28) 



where q is the total hadronic four-momentum. The spectral functions (amplitude form 
factors) are defined as sums of the corresponding Breit-Wigner amplitudes 



Fk{q) = Pi + PkBWp{q) + P',BWp,{q) + &Wp„{q) 



(29) 



the coefficients are functions of four-momenta describing the shapes of the corresponding 
contributions, the coefficients are complex amplitudes, the coefficients and (31 

are complex amplitudes of the fcth spectral function component, and the sum in Eqn. ( pSD 
is taken over all non-a; contributions, as described in Section [VII E| below. For practical 
purposes, the amplitude of one of the ak is set to a non-zero constant and the phase of one 
of the Ofc is set to zero; the corresponding rescaling factors are then included into the other 
amplitudes, while the overall normalization of the signal p.d.f. is set to unity. The explicit 
form of the form factor fj^ is based on Refs. which employs correct uj helicity 

factors and proper symmetrization. 

From the fit, we determine the amplitudes at for the decay r — > 27r^7r+7r°z/T-. We integrate 
these amplitudes over a large sample of r ^ Stttt^i^,- events to translate them into relative 
branching fractions. 



D. Background p.d.f. 

The total background in the selected data sample is estimated to be 4.9%. Because it is 
small, we do not calculate a p.d.f. for each individual component separately. Instead, we use 
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a simpler, empirical approach: we reweight the Stttt^ phase space to create an effective matrix 
element describing the shape of the sum of the r and qq backgrounds that was determined 
using detailed Monte Carlo simulations. We use a purely empirical form of this effective 
matrix element: 

Fb = {Ml - M(37r7r°)2) x (1 + 
0.8 0.8 



(M(7r%-)2 - M2)2 + r2 (M(7r07r2-)2 - M2)2 + r2 
0.011 0.011 ^ 



(M(7r07r+7rr)2 - M2)2 + Tl {MiuH+^^^f - M^f + Tl 
0.005 0.005 

TTTTZ^nZv) ?i,r2 , -n2 / 



The choice of the parameterization is largely ad hoc, the coefficients are tuned to describe the 
background distributions in seven mass projections, as shown in Fig. |^. These distributions 
are obtained from Monte Carlo simulations of r backgrounds and an appropriately scaled 
contribution from e^e~ — > qq continuum events. 



E. Fits to various models 

We use 4 models to parameterize the hadronic current. The models account for the 
following contributions: 

• Model 1: ujtt, fnm and non- resonant Stttt^; 

• Model 2: ujti and aivr; 

• Model 3: ujti, ain, crp and /o(980)p; 

• Model 4: utt, aivr and pTrn. 

Hadronic current expressions for each model are given in the Appendix. For the pim 
contribution, the spectral functions for the p^tt^tt^, p^Tx^ix" and p+Tr^vr^ components are 
assumed to be identical, and the corresponding form factors /f are parameterized in the 
same way as it is done in the TAUOLA package pO[. We use the PDG [0] value for the 
mass of the Oi resonance (Ma^ = 1230 MeV) and adopt the prescription of Ref. to 
use the higher end of the width interval recommended by PDG (Fq^ = 600 MeV). For 
the mass and width of o", i.e., /o(400 — 1200) in PDG's notation, we assume values ||2^ of 
M„ = 860 MeV and F^ = 880 MeV. We use PDG values for the mass and width of the /o(980) 
of Mfg = 980 MeV and F/^ = 70 MeV. To obtain the relative contributions of the r — *• STrir^Ur 
decay components, the fit parameters are integrated over 300,000 of Stctt^ phase space events 
and renormalized. Fit errors are propagated in accordance with this procedure. We use the 
following definition: Ri = Bi/B{T^ where Bi is a branching fraction for the 

corresponding channel, e.g., i?ai7r = B{t^ (ai7r)~z/T-)/i3(r~ — > 27r~7r"'"7r°z/.r)- Statistical 
errors on the relative fractions Ri are obtained from the maximum likelihood fits assuming 
one standard deviation. 
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FIG. 7. Background parameterization plotted for different submass projections. The predic- 
tions from Monte Carlo simulations are shown with crosses. The background parameterization is 
shown as a histogram. 
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To simulate the form factors BWpi, BWa^ and BW^r, we used both fits with mass- 
dependent and mass-independent widths. The observed variation in the amphtude values 
is small and has been included in the total systematic error. We fix the values of the mass 
and central width of p' in the fit; if these parameters are allowed to vary, the fit becomes 
unstable. 

To estimate the goodness-of-fit for each model, we generate 40 Monte Carlo samples 
with the matrix element averaged over the neutrino direction. Each sample is generated 
with parameters obtained in the fit to the data and the number of events in each sample 
is approximately equal to that observed in the data. The goodness-of-fit is estimated as a 
fraction of samples where —2Ln£ exceeds that in the data. 

Integrated amplitudes with statistical and systematic errors, as well as the confidence 
levels of the fits for the four chosen models, are listed in Table |V[ The calculation of the 
systematic errors is explained in the following section. Projections of the fits on various mass 
combinations are shown in Figs. |§-|TD|. 

Models 3 and 4 provide the best description of the data. Both these models are dominated 
by cuvr and Oivr with small additional contributions of the crp, /op or non-resonant pim 
channels. Submass projections for Models 3 and 4 are almost identical, and we show only 
those for Model 3. 



TABLE IV. 

Fit results for various models. Ri is defined as: Ri = Bi/B{T^ 27r^7r+7r°z/^), e.g., R, 
B{t- {aiTT)-v^)/B{T- 27r-7r+7rV^). 



Model 


Integrated amplitudes 


Sum of amplitudes 


Goodness-of-fit 


Model 1 


Rpo^-^o =0.11 ±0.01 ±0.02 
Rp-^^^+ = 0.19 ± 0.02 ± 0.04 
Rp+^-^- = 0.23 ± 0.02 ± 0.04 
R^^ = 0.40 ± 0.02 ± 0.05 
i?3^^o < 0.06 at 95% CL 


0.93 ± 0.04 ± 0.08 


15% 


Model 2 


R^^ = 0.38 ± 0.02 ± 0.02 
R^^^ = 0.43 ± 0.02 ± 0.02 


0.81 ±0.03 ±0.03 


< 5% 


Model 3 


R^^ = 0.38 ± 0.02 ± 0.01 
i^aiTT = 0.49 ± 0.02 ± 0.02 
Rap = 0.01 ± 0.02 ± 0.01 
i?/„p = 0.01 ±0.01 ±0.01 


0.89 ± 0.03 ± 0.03 


20% 


Model 4 


Ru;n = 0.39 ± 0.02 ± 0.01 
Ra,n = 0.50 ±0.03 ±0.01 
Rpo^-^o = 0.01 ± 0.01 ± 0.01 
Rp-^^^+ = 0.02 ± 0.03 ± 0.01 
Rp+^-^- =0.01 ±0.01 ±0.01 


0.93 ± 0.04 ± 0.02 


20% 



The presence of a non-resonant Stttt^ contribution does not improve the goodness-of-fit 
in the ci;7r-p7r7r-non-resonant-37r7r° model, and we set an upper limit on the non-resonant 
contribution shown in Table 0. We also ignore the non-resonant 37r7r'' component in the 
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two other models. 



F. Systematic errors 



To estimate systematic biases and errors for every model, we generate 40 Monte Carlo 
samples with amplitudes obtained from the fit to the data and then fit them to the chosen 
parameterization of the hadronic current. Each sample is mixed with simulated background 
events in the proportion expected from the data. We integrate the results of each simulation 
to estimate necessary corrections and systematic errors for each of the amplitudes. Biases 
of the reconstructed amplitude values are shown in Table |IV|. Errors are determined as root 
mean squares of the parameter fluctutations. 

The total systematic error also includes uncertainty due to the choice of the parameters 



of p' . We use PDG [15| values for the masses and widths of the p and p" mesons and vary the 



p' parameters within the errors obtained in the fit described in Section |rV|. We also switch 



to a smaller mass of p' |jTT|: Mpi = 1370 MeV/c . The results appear to have very little 



sensitivity to the choice of the parameters for p'. 

All contributions to the systematic error are added in quadrature. Systematic biases for 
the components of the hadronic current are shown in Table 0. 

TABLE V. 

Systematic biases for the hadronic current components. Systematic bias is defined as 6i = 

^-p^reconstructed _ j^generated-^ j j^generated 



Model 


Systematic bias 6i 




6pO^-^o = (+4.0 ± 3.8)% 




6p-^-^+ = (-7.6 ±3.8)% 


Model 1 


6p+^-^- = (+2.3 ±2.6)% 




= (+2.5 + 2.6)% 




S^^^o = (+1.2 + 9.0)% 


Model 2 


= ( + 1.0 + 2.5)% 




Sa,^ = (-2.0 + 2.5)% 




6^^ = (+1.0 + 2.5)% 


Model 3 


Sa,n = (-2.0 + 2.6)% 




6ap = (0.0 + 14.0)% 




^foP = (+1-0 ± 16.0)% 




Su.n = ( + 1.0 + 2.5)% 




5a,^ = (-2.0 + 2.5)% 


Model 4 


Spo^^^o = (1.0 + 23.0)% 




= (-0.0 ±15.0)% 




6p+^-^- = (-2.5 + 33.0)% 
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FIG. 8. Comparison between fits and data distributions for the wtt, pir-K and non-resonant 
Stttt^ model. The horizontal axis is given in units of GeVjc? . Solid squares represent the data, 
solid line represents the nominal fit, and the dashed line represents the parameterized background. 
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FIG. 9. Comparison between fits and data distributions for the lott and oitt model. Tlie 
horizontal axis is given in units of GeV/(?. Solid squares represent the data, solid line represents 
the nominal fit, and the dashed line represents the parameterized background. 
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FIG. 10. Comparison between fits and data distributions for the lott, oitt, ap and fop model. 
The horizontal axis is given in units of GeV/ (? . Solid squares represent the data, solid line represents 
the nominal fit, and the dashed line represents the parameterized background. 
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G. Comparison with similar analyses 



A recent analysis by the ALEPH Collaboration [Q assumes ujtt and non-resonant pmr 
contributions and uses a simplified fit to the submass projections without assigning sys- 
tematic errors. This analysis by ALEPH gives the following amplitudes: RpO^-^o = 0.11, 
Rp-^-^+ = 0.20, Rp+^-T,- = 0.22, and R^t^ = 0.40. Predictions from the chiral perturbation 
theory p3[ for the same model give: Rp^^^^+ = 0, Rp+^-^- / RpO^-^o = 2. These results can 
be also compared with the "standard mix" expectations: i?pO^-^o = 0.28, Rp-^~^+ = 0.27, 
Rp+^-^- = 0.08, and i?^^ = 0.30. 



H. CVC test 

The CVC hypothesis relates the weak charged hadronic current in the r lepton decay 
to the isovector part of the electromagnetic current. Thus, the spectral function (H) can be 
related iTBIl to the cross-section for e^e~ — > 47r as: 



V'-\q) = ^,<ye.e-^U<l) . (31) 

Therefore, the resonant structure of the hadronic final state in the decay r 37r7r°z/^ should 
be comparable to that in the reactions e^e" 27r+27r^ and e^e' 7r+7r^27r'^ in the energy 
range corresponding to the r mass. The only expected difference is due to the fact that the 
cuTT final state is not accessible in the all charged final state. A recent analysis of the CMD-2 



data |25| shows clear dominance of the oivr component in the four charged pion final state 
and of the aivr together with the ujn components in the 2 charged and 2 neutral pion mode. 
A plot comparing the 37r7r'' and ujtt spectral functions in the decay r 37r7r*'z/^ with the 



corresponding cross-sections measured by the CMD-2 group is shown in Fig. |1T| Our 
results are in excellent qualitative agreement and support the CVC hypothesis. 



VIII. CONCLUSIONS 



The large data sample of r ^ 37r7r°z/^ decays collected by the CLEO II experiment 
allowed for a detailed study of the resonance decomposition and structure functions of the 
37r7r° final state. The analysis indicates dominance of the uott and aivr contributions which is 
in good agreement with theoretical expectations based on the CVC hypothesis and Chiral 
Langrangian calculations. At the same time, with the present statistics, other models can 
provide satisfactory descriptions of the existing resonant structure. The goodness-of-fit is 
slightly higher for the models with a dominating oivr contribution, though the Oivr channel 
alone does not reproduce the observed non-uj contribution precisely. The mass and width of 
the p' resonance have been extracted from the fit to the r — >■ cuttz/^ spectral function, and a 
new upper limit on the non-vector current contribution to this spectral function of 5.4% at 
90% CL has been obtained. 



30 



0.16 



^ 0.12 

CM 
CM 

^ 0.08 



0.04 - 



T 1 r 



0.90 



T 1 r 



T 1 r 



T 1 r 



□ ■ ^ A ^ 

g 5 A 



0.00 ^ ^ L 



J L 



1 



A A 



A' 



J L 



1 



J L 



1.10 1.30 1.50 

M(47t), GeV/c' 



1.70 



FIG. 11. Comparison of V'^'^'^" (q) and V'^'^{q) measured in this analysis with the corresponding 
cross-sections measured by CMD-2. Squares represent the full four-pion spectrum, and triangles 
represent the utt component. Our results are shown with open symbols, and the results obtained 
by CMD-2 are shown with filled symbols. 
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APPENDIX A: HADRONIC CURRENT FOR THE DECAYS r 

AND T Fopur 



Expressions for the hadronic current for r decays into the four pion final state via aivr, 
ap and /op channels are derived under the assumption that the r decays into an interme- 
diate vector resonance p: r pz/^, which consecutively decays into the shown channels. 
The vector boson p has quantum numbers I{J^^) = 1(1 ), and its spectral function, in 
accordance with Eqn. (^), is simulated as mixture of p and p' resonances and a phase space 
contribution. Assuming isospin conservation and enforcing Bose symmetry, one obtains the 
following expressions for the form factors /f in the decay — > 27r^7r+7r°z/T-: 

n,{Qi)BWa,iQl) [BW.iiq, + g2)')(gi - ^2)" + BWpHq, + q^f){q^ - q^Y 
- TUQ2)BWaAQl) \BW,{{q, + q,f){q, - q,Y + BWp{{q2 + q,f){q2 - q,Y 



- TUQz)BWaAQl) BWp{{q, + g3)')(gi - q^T + BW,{{q, + q,f ){q, - q^Y 



(Al) 



fl^=T^^{q)^ 



where 



BW,{{q^ + q2f)BWp{{q^ + q,fYq, - q^Y 

+BW„{{q, + q4Y)BW,{{q2 + q^YYq2 - qz) 



(A2) 



qi 
q2 
q-i 

?4 



four-momentum of vr"*" ; 
four-momentum of 7rf ; 
four-momentum of 7r° ; 
four-momentum of vr^ ; 



(A3) 



g = gi + 92 + gs + <?4 ; 
Qi = qi + q2 + q4 ; 

Q2 = gi + g2 + gs ; (A4) 
Qs = gi + gs + g4 ; 

and T^"{x) denotes the projection operator: 

T'^-(x) = -g^^- + —— . (A5) 

The expression for the r decay into four pions via the /op channel is identical to Eqn. ( [A^ ) 
with the mass of /o substituted instead of the mass of a. 
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